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Abstract 


A higher order theory which satisfies zero transverse shear stress conditions on 
the bounding-planes of a generally laminated fibre-reinforced composite plate subjected to 
transverse load is used to incorporate with the Von Karman’s geometric nonlinearity. A C° 
continuous displacement finite element formulation with seven degrees of freedom, is pre- 
sented and the coupled membrane-flexure behavior is investigated. The 9-noded Lagrangian 
elements are used for mesh elements of the plate.The results for isotropic, orthotropic and 
cross-ply plates are found. The material non-linearity is also introduced to make a complete 
nonlinear study for bending of the laminated plates.The variation in bending behavior is ob- 
tained with a/h by using linear HSDT, geometric nonlinear HSDT and completely nonlinear 
HSDT. Some cases of angle-ply laminates are also presented. The first-ply failure analysis for 
various lamination-schemes is also done for transverse loads only, and results are compared 
with those available in various literature. 
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Chapter 1 


Introduction and Literature Review 


I 


1.1 General 

In recent years, advanced composite materials are being used increasingly in many en- 
gineering and civilian applications, ranging from the fuselage of an aeroplane to the frame of 
tennis racket. The high stiffness-to-weight ratio, coupled Avith the flexibility in the selection 
of lamination scheme that can be tailored to match the design requirements, make a lami- 
nated plates attractive structural components for automotive and aerospace vehicles. The 
increased used of laminated plates in various structures has created considerable interest in 
their analysis. 

The analysis of plate and shell structures is of considerable interest in various areas 
of structural mechanics. It is, therefore, natural that with the development of the finite 
element method a large number of different finite elements have been formulated for the 
analysis of plate and shell problems. In these developments, basically two approaches have 
been followed: firstly, what maybe called a classical procedure and, secondly, an approach 
in which displacement/rotation isoparametric elements are employed. 
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1.2 Literature review 

The two-dimensional analyses of laminated composite plates in the past have been based 
on one of the following two theories: 

(1) the classical lamination theory, 

(2) shear deformation theories. 

In both these theories it is assumed that the laminate is in the state of plane stress, the 
individual lamina are linearly elastic, and there is perfect bonding between layers.The clas- 
sical laminate theory, which is an extension of the classical plate theory(CPT) to laminated 
plates, ignores the transverse stress components and models a laminate as an equivalent sin- 
gle layer.The first complete laminated anisotropic plate theory is attributed to Reissner and 
Stavsky[l] .The classical laminate theory is adequate for many engineering problems.However, 
laminated plates made of advanced filamentary composite materials, such as graphite-epoxy 
(whose elastic modulus to shear modulus ratios are very large), are susceptible to thick- 
ness effects because their transverse shear modulii are significantly smaller than the effective 
elastic modulus along the fibre direction.These high ratios of elastic modulus to shear mod- 
ulus render the classical laminate theory inadequate for the analysis of composite plates.An 
adequate theory must account for accurate distribution of transverse shear stresses. 

To improve the situation, many theories which account for the transverse shear and 
normal stresses are available in the literature. The Reissner[2]-Mindlin[3]-type theories for 
isotropic plates was extended to laminated structures [4] .This theory is normally called the 
First Order Shear Deformation Theory or FSDT. These theories are based on an assumed 
displacement field that is expanded in terms of the thickness coordinate.A generalisation 
of the FSDT for homogeneous isotropic plates is due to Yang et al. [5] and Whitney and 
Pagano[6].In this theory, the normals to the mid-plane before deformation remain straight 
but not necessarily normal to the mid-plane after deformation,and consequently, a correction 
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to transverse stiffnesses is required. Since the transverse shear strains are constant through 
the thickness, the transverse stresses are also constant. And therefore, these shear deforma- 
tion theories do not satisfy the condition of zero transverse shear stresses on the top and 
bottom surfaces of the plate, arid require the shear correction to the transverse shear stiff- 
nesses. Introducing the shear correction factors, whose accurate prediction for anisotropic 

laminates is cumbersome and problem dependent, J.M.Whitney[7], is not always desirous 

i 

and pleasing. 

Thus in order to have a reliable analysis and safe design, the proposal and developments 
of models using higher order shear deformation theories have been considered. Lo et al.[8] 
and [9] reviewed the available work in the field and formulated a theory which accounts for 
the effects of transverse shear deformation, transverse strain and polynomial distribution of 
the in-plane displacements with respect to the thickness co-ordinate. The theory showed 
that for problems which involves rapidly fluctuating loads with characteristic length of the 
order of the thickness of plate, one of those polynomial theories is required to obtain the 
meaningful results. 

Third order theories have been developed by many theoreticians e.g. Reddy [10], Li- 
brescu[ll], Schmidt[12], Murty[13], Levinson[14], Seide[15], Murthy[16], Bhimaraddi and 
Stevens[17], Mallikarjuna and Kant[18], Kant and Pandya[19], and Phan and Reddy[20] 
among many others. Third order theories are simplest of the all the higher order theories, 
yet the desired degree of accuracy is very good.So the third order theories are very popular 
when computational complexities are desired to be avoided. In higher order theories, an 
additional dependent unknown is introduced into the theory with each additional power of 
the thickness coordinate. In addition,these shear deformation theories do not satisfy the 
condition of zero transverse shear stresses on the top and bottom surfaces of the plate, and 
require a shear correction to the transverse shear stiffnesses. The three dimensional theories 
of laminates, in which each layer is treated as homogeneous anisotropic medium(see Ref. [21]) 
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are intractable as number of layers become large. Thus, a simple two-dimensional theory 
of plates that accurately describes the global behaviour of laminated plates seems to be a 
compromise between accuracy and ease of analysis. 

The displacement field used by Reddy[22] is similar to that of Levinson[14] and Murthy[[16], 
but the latter uses the equilibrium equations of the classical plate theory, which are inconsis- 
tent with the assumed displacement field.Reddy’s modifications consist of a more systematic 
derivation of displacement field and variationally consistent derivation of the equilibrium 
equations.The conventional variational formulation of the classical plate theory as well as 
higher-order theory involves higher order(i.e. second order) derivatives of the transverse 
displacement.Therefore in finite element modelling of such theories one should impose the 
continuity of not only the transverse displacement but also its derivatives along the element 
boundary.In other words, a conforming plate bending element based on displacement formu- 
lation of these theories requires continuity of transverse displacement and their derivatives 
across the inter-element boundaries [23] .The construction of such an element is algebraically 
complicated requiring, for example, a quintic polynomial with 21 degrees of freedom for a six- 
node triangular element .Computationally the element requires much storage and computer 
time. 

To overcome the stringent continuity requirements placed by the conventional varia- 
tional formulation of the classical plate theory , several alternative formulations and asso- 
ciated elements have been proposed (see (Putcha and Reddy[24])). These include the hybrid 
fini te elements and mixed finite elements.The hybrid elements are based on variational state- 
ments that use independent variation of displacements inside the domain of the element and 
tractions on the boundary of the element.The mixed elements use stationary variational 
principles, such as the Reissner variational principle or Hu-Washizu variational principle, to 
construct independent variations of both displacements and bending moments in a plate.The 
conventional variational formulation of the FSDT leads to a C°-element, often referred to as 
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the Mindlin element.Thus Mindlin. plate elements are particularly easy to formulate. Such 
formulations readily accomodate moderately thick or thin plates of variable thickness, curved 
boundaries and laminated plates.The initial confidence expressed in these elements was later 
subject to doubts on two accounts.First, the serendipity elements were found to exhibit an 
overstiff locking behaviour in very thin plate situations and secondly, the Lagrangian ele- 
ments with reduced integration rules contained spurious zero energy modes not associated 

i 

with rigid body movement.However, extensive studies (Pugh et ol.,1978) on this phenomenon 
led to the conclusion that the nine-noded Lagrangian element with selective integration of 
the shear terms is the best for such elements. 

The present study uses the Higher Order Shear Deformation theory of Kant and Pandya[25] 
But this study differs from the previous theories in two ways: 

(1) it incorporates Von Karmans’ non-linearities terms,and 

(2) it considers non-linear constitutive laws too. 

The geometric non-linearity is incorporated into the conventional HSDT formulation pro- 
posed by Kant et a!. [25] .This makes the finite element formulation nonlinear. The bound- 
ary conditions are of mixed type and the governing differential equations are non-linear.A 
nine-noded Lagrangian C^-continous element with seven degrees of freedom per node is pro- 
posed. The solutions to problems with material nonlinearity were attempted by numerous 
researchers (see Ref. [26]), for shells and plates.The analysis of laminated composite plates 
with material nonlinearity was done by Reddy et a/. [27]. In which, the modified Ramberg- 
Osgood[28] formulation is utilised for the calculation of tangent modulus.In another work, 
a cubical formula is proposed for the shear stress calculation only, and linear stress-strain 
relation is used for the normal stress calculations[29]. 

Petit and Waddoups[30] simulated the nonlinear stress-strain response of a lamina by 
using an increment method. Hahn and Tsai[31], by using complementary elastic energy 
density, presented a constitutive equation for in-plane shear stress and strain of a composite 
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lamina. Sandhu[32], by using cubic spline interpolation functions, represented experimental 
shear stress-strain data. Jones and Nelson[33] modeled nonlinear mechanical properties as 
functions of strain energy density. Jones and Morgan[34] showed that Jones-Nelson model 
converged only up to a specific strain energy values.' 

In the present work, we utilized both of these approaches [28] and [29], so that, for the 
normal stresses,the modified Ramberg-Osgood formula is used and for taugential stresses, 
the cubical stress-strain formula is used upto the limit of practicality of stresses. Beyond the 
limit of practicality of the stresses, the structure is seldom encountered with larger stresses. 



Chapter 2 

I 

i 

Model Problem 


2.1 Kinematics: 

The development of the present HSDT begins with the assumption of the displacement 
field in the followin^g form: 

i 

u{x, y, z) = Uo{x, y) + z* 9x{x,y) + z^ * ul{x, y) + * dl{x, y) (2.1) 

v{x, y, z) = Vo{x, y) + z* Oyix, y) + z^* Vq{x, y) + z^* 9*{x, y) (2.2) 

w{x,y,z) = WQ{x,y) (2.3) 

where uq, vq and wq denote the displacements of a point (x,y) on the mid-plane and 6{x) 
and 9{y) are the rotations of the normals to the mid-plane about y- and x-axes respectively. 
The parameters uI,Vq, 9% and 9* are the corresponding higher order terms in the Taylor’s se- 
ries expansion and are also defined at the mid-plane. The condition that the transverse shear 
stresses vanish on the top and bottom faces of the laminate is equivalent to the requirement 
that the corresponding strains be zero on these faces. 

The transverse shear strains are given by: 
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^ , 9(wl a + 0 + ^ _l o + -v2 * a* , 9(ii>o) 

(2.4) 


i 

i 

Jxz—q(^2 ) ■*■ a(x) + 2 *z*Uq + S*z * 6j.+ 

(2.5) 

Now, equating, 

7yz =7xz = 0 at(x, y, ±{h/2)) 


we get. 

\ 

Vq = uI = 0 

(2.6) 

and 



3/^2 + Qi^y) ) 

aw ' 

using these relationships, we get the new displacement field as: 

w = uo + ^ 

(2.7) 

■ (2.8) 

(2.9) 


V = Vo + z * 9y + * 6* 

(2.10) 
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w = wo (2.11) 

Now, we shall proceed with displacement field given by equations(9,10,ll ), and later 
we shall introduce conditions (7,8) into the shear-rigidity matrix,!)^. The strains associated 
with the displacement field are: 


^xq ^ ^xx “i* - 2 ? ^ ^xx 

(2.12) 

“ ^yo + -2^ * ^yy + -2^ * Kyy 

(2.13) 

= Jxyo +Z*Ka:y + Z^^ K^y 

(2.14) 

lyz ~ ^yy + - 2 ? * (f>yy 

(2.15) 

'Jxz ” ^xx “h * (l>xx 

(2.16) 


Here, we shall introduce the Von Karmans’ non-linearity terms to formulate the geometric 
non-linearity problem.Therefore, 


dun 1 ,dw.n T 


NL 

’XX 


_ ^ 1 i ^ NL 

dy^ 2 w + Sy 

_ duo ,dvo ^ 

'7x1/0 - dy dx^ dx dy 




(f^x 


Ox + 


(■* 

"xx 


dwp 
dx ' 


dOx^ 


dOy 


_ dOx 


dx ’ 

^yy 

dy ’ 

1 ^xy 

dy 


d0*x. 

= 


: kZ.. 



dx’ 

^yy 

dy' 

' *''xy 

dy 

dx 

<l>yy = 

■ Oy + 

dwo 

dy’ 

^Ix 

= 3 . 0 ^ 



(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 
( 2 . 22 ) 
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Now, the strain vector can be written as, 


f ^xx 1 


^yy 


{^}“‘ I Jxy f 


; 7yz 

1 , 


{e^} + {e^^} + 2 : * {«} 4- * {^*} + * {k*} + {<f>} 


(2.23) 


where, 


{£‘}={ 


^xo 5 ^yo ) Txvo ’ 


•,1 

0, 0} 


{f“}={ 0.5 . (t)^ 0,5 . (^)^ (t * t). 0, 0 }’ 


T ^ ^ 

■[/c}=| Kxxi l^yyi ^xyi 0) 0 ^ i } | ^yyi ^xy’ ® J 


(2.24) 


(2.25) 


(2.26) 


{^}— I 0, 0, 0, <l>yyj ^xx } ) ^yy’ } 


(2.27) 
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2.1.1 Constitutive equations: 

The constitutive equation for the layer can be written with reference to the lamina 
co-ordinates (1,2,3) as 


> 

0-1 

L 

' Cn 

Ci2 

0 

L 

r N 

<72 

> = 

Cl2 

C 22 

0 

* ^ 

C2 ' 

Ti2 

\ J 


0 

0 

1 

CO 


7i2 

V. y 



(2.28) 


(2.29) 


The following relations hold between these and the engineering elastic constants: 


Cii = — 

1 - 1/12 * 1^21 

(2.30) 

1/12 * El 

012 ~ 

1 ~ *^12 * ^^21 

(2.31) 

Eo 

C 22 = 

1 ~ *^12 * ^^21 

(2.32) 

C33 = Gi 2 

(2.33) 

CiA = G 23 

(2.34) 

C'55 = Gi 3 

(2.35) 
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Now, in the laminate co-ordinates (x,y,z), the stress and strain relationships for the 
lamina are: 


|' ^xx 

L 

Qll Qi 2 Qi 3 

L 

f \ 

^xx 

j 

► = 

.Qi 2 Q22 Q23 

* < 

^yy 

i . 


Qi 3 Q23 Q33 


Ixy 

V, ^ J 


(2.36) 



Qu 

Q45 



where, 


W={ o-j: 


'yyi *yzj 


(2.37) 


(2.38) 


{^} { ^yyi '1xy> 'lyzi Ixz } 

(2.39) 

are stress and strain vectors with respect to laminate axes and Q^ -are the plane stress 
reduced elastic constants in the laminate axes of the V'^ lamina. 




2,1 Kinematics: 


13 


2.1.2 Material model: 

If the material used in the laminate construction is non-linear, i.e. the stress and 
strain relationships are non-linear, then the co-efficients Qij vaxy with the deformation.If it is 
assumed that the nonlinearity can be represented by uncoupled stress - strain behaviour for 
each of the strain components, then a suitable curve fitted to the experimentally observed 
stress - strain behaviour can serye as the material constitutive law. The Ramberg - Osgood 

I 

relationship is one such procedure. 

For each of the normal strain components , the stress - strain relationship can be simu- 


lated by the Ramberg - Osgood relation as: 

a . 




(2.40) 


where e is the strain, a is the stress, jFf is the initial modulus, and A and m are the 
parameters defining the curve as derived at two different values of strains in the stress-strain 
curve. ! 

A modified stress - strain relationship can be derived from the above equation (2.40) as 
in ([27]), 

(2.41) 




where Ep is the elasticity modulus at the plastic range and Op is the stress at which the 
tangent to stress-strain curve, at some value of strain in the plastic range, hits the stress- 
axis. Ep is the angle of this tangent with strain-axis. Tangent modulus can be defined 


Et = ^ = 


{Ei - Ep) 


[ 1 +^ 


I (m+l)/m 


(2.42) 


For the tangential strain components, the stress - strain relationship can be written as, 
for the smaller values of strains i.e. upto around 2.5 percent,: 

Txy — Gxy * 'Ixy Afi ♦ ^xy Oij/ 


(2.43) 


2.2 Total potential energy and variational formulation: 
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Tyz Gyg ♦ ')yg A 4 sl" Gyz * 


yz 


'^xz Gxz * ')/xz <^5 * Gxz * T; 


,3 

xz 


(2.44) 

(2.45) 


in which Ae, A 5 and A 4 denote the degree of nonlinearity with different shear strains and can 
be found froni experimental curves [31] to suit the purpose of applicability. 

The tangent moduli for these sheair curves are: 


where, 



= (/- 

-3* A' 

*g' * 7 ^* 

(2.46) 

{^t) = { - 



Et.. I"" 

(2.47) 

{t'} = -1 

L '’’^1 


r..f 

(2.48) 

{V} =i{ 

Ixyt 

lyz’i 

lx, } 

(2.49) 

{<?"}=( 

Gxyj 

^yz} 

Gxz } 

(2.50) 

{A'} = 

{ 

^5? 

A 4 1 

(2.61) 


Equations (2.42) and (2.46) are used for each strain level to update the stifiEness matrix, 
which shall be found in finite element formulation. 


2.2 Total potential energy and variational formula- 
tion: 

The first variation in total potential energy 11 of the given system should be zero for 
the equilibrium of the system: 


(jn = o 


(2.52) 
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The displacement vector d is defined as: 


3. 

{d}={ uo, uo, Wo, 0*, Oy $1, } 


(2.53) 


The functional given by total potential equation is minimized while carrying out the explicit 
integration through the thickness.This leads to the following variational formulation: 

dV- Jj(f • P dA = 0 (2.54) 


The internal virtual work is: 


SUi 


i = / Se^ • a 
Jv 


dV 


or, 


* ■^x * ^^y ^^xy * ^^xy ^^xx * ^^x ^^yy * ^^y 

SKly . M*, + 5<f>.. • Q, + 6<l>yy . Qy + 5^, • Ql + 5cl>ly • Q;] dA 

where the following 13 stress - resultants for the n-layered laminate are given as: 

dz 



^x 

II 

Cfx 

< 

Ny 



Nxy 

j 

^ '^xy 

' : 

Ml ' 

-^Cr 

CTx 

< My : 

m; 





Txy 

f Qx : 

Ql ] 

q; ) 

-It. { 

'^XZ 1 

V 


[2 : 7^^- dz 


(2.55) 


(2.56) 


(2.57) 


(2.58) 


(2.59) 
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After the integration of these equations and using stress - strain relationships, 
matrix notation: 

f JV 1 - 

M IB] [0] 


\ 

K 

P 

M 


Ml 

► = 

Q 


Q* 

j 



Kl [A1 [0] 


1 

^ V 

€ 


K 

< 




-* 

J’ . 


> • dz 


in which, 


JV = { JV,. Ny, 

m = m„ f ;Ar = { m;, m;, 

<?={ <3., <3, r;«'={ «, e: Y 

^ ^ ^xxi 'Yxy } { l^xxi ^yyj ^xy } ^yy'> ^xy J 

^ { 4 ^xxi ^yy } > 4^yy } 


and. 


[■^laxa ~ 

L=1 


QiiHi QnHi QizHi 

SYM. QnH^ Q 23 H 1 

QzsHi 


£th layer 


[■^Isxe ~ X) 

L=1 


Q11H2 Q12H2 Q13H2 Qi\H^ QnH^ Q1ZH4, 

Q12H2 Q22H2 Q23H2 Qi2Hi Q22-H4 Q23-H4 

Q13H2 Q23H2 Q33H2 Q13H4 QtzH^ Q33H4 


£th layer 


we have in 

(2.60) 

(2.61) 

(2.62) 

(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 


(2.67) 
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PJ6.6 = E 

L=1 


QiiHz 

Q 12 H 3 

QxzHz 

QiiHz 

Q 12 H 5 

Q 13 H 5 


QuHz 

Q 22 HZ 

Q 23 H 3 

Q 12 H 5 

Q22H5 

Q 23 H 5 


Q13H3 

Q23H3 

Q33H3 

Q13H5 

Q23H5 

Q33H5 


QnHz 

Q 12 H 5 

QxzHz 

QxxH^ 

Qx2Hr 

QxzHj 


Qx 2 Hz 

Q 22 HZ 

Q 23 H 5 

Qx2H7 

Q 22 H 7 

QizHj 


QxzHz 

Q 23 HZ 

QzzHz 

QxzHt 

Q23Hr 

Q 33 H 7 


jjkh layer 


[AI4.4 = E 

L=X 


Q 55 H Q^zH 0 0 

Q^aH 0 0 

Sym. QzzH* Q^zH* 

QuH* 


jth layer 


wherein; 


\{h\ — hi_i) 


= h 

— ^l-x) 

= iHi-Hzx4/h‘^) 


( 2 . 68 ) 


(2.69) 


(2.70) 


Hx = (hL-hL-x) 

H 2 = m-hi-x) 

Hz 4 m-hl-x) 

Ha = m-hU) 

Hz 
H 7 
H 

H* = (Hz-HzxhyA) 

where h is the overall thickness and hi is the thickness of Lth layer. 

The shear-rigidity matrix is evolW by incorporating an alternate form of the equations 
(2.15) and (2. 16), viz., 

^yy \^yy ~ ^ 

(2.71) 

<i>xx + ~ ^ 


and using these relations (2.71), the resulting theory becomes consistent in the sense that it 
satisfies zero transverse shear stress conditions on the bounding-planes of the plate. 
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2.3 Governing equations cind boundary conditions: 

Prom equations (2.24)-(2.27) and (2.56), we have 

= ^[{(<^^0,* • Nx + SUQ^y • Nxy) + ((5U0„ • Ny + (5U0.. • N,;y)} 

+{{Nx • 'Wo.. • 5wo,^ + Ny • wo,« • Swq^^) 

+{Nxy ■ 'Wo,. • (Jwo,„ + Nxy • 5wo^ + Qx • ^wo,., 

+Qv • ^Wo,y)} + {{Mx • + Mxy • 501, + Qx • <50x) 

+ (Mj, • 50J,,, + Mxy • 50y,,, + Qy ' 59y)} (2.72) 

+{(K* • sei. + M:y • 50;, + 3 • q: • 50:) 

+(m; . 50;, + . 50;, + 3 • Qise;)}] dA 

And , the external virtual work is: 

SWe=fS(f-PdA (2.73) 

J A j 

In equation (2.73) , the load vector P consists of two parts, viz., 

(i) transverse load, P,, and 

(ii) in-iplane loads, Px, Py, Pxy 

In bending problem , we shall consider only the transverse loads. Therefore , load- vector 
is: 


{P}={ 0 0 P, 0 0 0 0 } 


T 


(2.74) 

and displacement vector is given by equation (2.53). Therefore , from equations (2.53), (2.73) 
and (2.74), we have 

5We= f Swo-PgdA 

JA 


(2.75) 
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Therefore , for equilibrium 

6Ui-5W,=^0 (2.76) 

+ Suo^y ■ Nxy) + (5uo,^ • Ny + (Juo,x • Nxy)} 

+{{Nx • Wo,, • 5wo„ + Ny • Wo,^ • 5wo,^) 

+iNxy • wo„ • + Nxy • wo„ • 5wo„ + Qx • Swq^^ 

+Qy • Swq^^)} + {{Mx • S9x^ + Mxy • S9xji + Qx • 

+ {My • 50,,^ + M:,J, ■ 59y,^ + Qy * 59y)} (2.77) 

+{(M: • . 50*,^ + 3 . Q: . 591) 

+(K • K + 3 • - iP<^ • Swo)} dA = 0 


Integrating the expressions in equation (2.77) by parts, and collecting the coefficients of 


Swo,69x,S9y,S9l and 69*, we obtain the following euilibrium equations in the domain 


n : 


Suq 

6vo 

6wo 


< 


69x 

69y 



ai\fx 

dx 

dMy 

dy 

dMj 

dx 

du; 

dy 


+ 

+ 

+ 

+ 



+ Qx 
+ Qy 


2Mh 

dy 

dM'„ 

dx 


-3-Q: 


3-q; 


+ Ny 



0 

0 


0 

0 

0 

0 

0 


(2.78) 


The boundary conditions shall be of the form: 
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Uo 

^0 


Wo 

Ox 

Oy 

0*x 


or 

or 


or 


or 


(W.+WJ 
(Ny + N^) 
Ny 

+Ny 


dw(i 
dx 

dy 

+^) 


\ "^Qx + Qy ) 

{Mx + ^xy) 


or “1“ iWjpy) 

or {Ml+Mly) 
or + 


( 2 . 79 ) 


This completes the derivation of governing equations.In the present investigation we 
have taken two types of boundary conditions, e.g. (i)simply-supported boimdary at all 
edges and (ii)clamped boundary at all edges.The boundary conditions as in figure 2.1 are 
shown only for the quarter plate due to symmetry of the plate. 
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2.4 Failure criteria 

The criteria for the failure analysis of laminates, at macro level, which is ply level failure, 
can be classified into two groups, namely (i) independent failure criteria, and (ii) polynomial 
failure criteria. In independent failure criteria the interaction among different stresses and 
strains does not exist. Therefore, in general, the polynomial failure criteria proposed by 
Tsai-Wu (1971) is used. The tensor equation of these failure criteria is expressed as: 


FijaiCj + FiCi > 1 (2.80) 

where Fi and are the strength tensors of second and fourth rank, and (jj (contracted 
stress notation is used) is the stress component in the material coordinate system. In a more 
explicit form, the criteria (2.80) is 


Fi(7i + F20’2 + 2Fi2(Ji<7’2 + 2Fiz<^l<^3 + 2F23(T2(T3 + FnO^ 

+F 220'2 + Fzzof + ^ 44(74 + Fsstrl + Fq^Oq + ...>! (2.81) 

in which ai, a 2 and crz are the normal stress components and in three principal material 
directions, while 0 - 4 , 0-5 and are the shear stress components at a point in planes 2-3, 1-3 
and 1-2 respectively. 

Following are the various polynomial criteria used (Reddy and Pandey, 1987 and Reddy 
and Reddy, 1992, (see Ref.[36]) for details): 

• Maximum Stress Criteria 

• Maximum Strain Criteria 

• Tsai-Hill Criteria 


• Hoffman Criteria 


5.5 Definition of failure 


23 


• Tsai-Wu Criteria 

Composite materials are made by mechanically very dissimilar phases; namely a stiff 

elastic brittle fibre and a compliant yielding matrix. The fibre strength is direction depen- 

-%■ ’ 

dent; and it possesses greater strength in the longitudinal than in the transverse direction. 
Hence, the strength of a lamina along the fibre direction is mainly due to the fibre strength 
and transverse to the fibre direction is mainly due to the matrix. Due to difference in inter-i 
laminar stresses, fibre debonding is also a possibility. Consequently, laminates fail in variety 
of different modes, but all failure modes can not be represented by single failure envelope. 
For convenience, the failure modes are categorised into the following groups: 

Fibre Mode: Tensile fibre mode is described by breakage and compressive fibre mode 
is described by fibre buckling. 

Matrix Mode; Matrix mode of failure is described by the failure due to in-plane 
normal stress transverse to the fibre direction, in-plane shear stresses and transverse shear; 
in the lamina. 

Delamination Mode: It is produced by transverse interlaminar stresses or transverse 
normal and shear components in layers adjacent to the interface. 

2.5 Definition of failure 

2.5.1 Lamina/Ply failure 

The ply failure is said to have occurred when the state of stress at any point within the 
lamina satisfies the selected failure criterion. The first-ply failure refers to the first instant 
at which one or more plies fail at the same load. 

2.5.2 Onset of delamination 


The following two criteria are employed to predict the onset of delamination: 


2.5 Definition of failure 
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• The maximum stress criterion 

• Interlaminar failure criterion([36]) 

As for the maximum stress criterion, delamination at any interface is said'%o have occurred 
when any of the transverse stress components in any of the two layers adjacent to the inter- 
face becomes equal to, or greater than its corresponding strength. As for the interlaminar 
failure criterion, the onset of delamination takes place when the interlaminar transverse 
stress components satisfy the following expressions: 

+ (2.82) 

where is the transverse normal stress component; ri 3 , r 23 are the transverse shear stress 
components in the principal material planes 1-3 and 2-3, respectively; ann is the peel 
strength and a^s is the interlaminar shear strength; these are taken equal to the tensile 
normal transverse strength and transverse shear strength ( corresponding to the plane 1-3) 
of a lamina, respectively. 


Chapter 3 

Finite Element Formulation 


3.1 Finite element formulation and solution technique: 


T he standard finite element technique is followed, in which the total solution domain 

1 

is descretized in ’NE’ sub-domains (or elements) such that 

NE 

n = £n‘{<i) (3.1) 

e=l 


where II and 11® are the total potential of the system and element respectively. The 
element potential can be expressed in terms of internal strain energy and the external 
work-done W® for an element ’e’ as: 


n®(d) = 17® - W® 


(3.2) 


in which d is the vector of unknown displacement variables in the problem and defined by 
equation (2.53). If the same interpolation functi on is used to define all the componenets of 
the generalised displacement vector d , we can write, 


nn 

d = Tij ♦ dj 
i=l 


(3.3) 


where , Ui is the interpolating (shape) function associated with node i ,di is the value of d 
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corresponding to the node i and nn is the no. of nodes in the element. 
Let us say the degrees of freedom are given as; 


Now, the variational energy is: 


tio = * n] 

(3.4) 

Vo = 53 ^^? * 

(3.5) 

Wo = * n? 

(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


i 


sn^id) = 5U‘ - dW^ (3.11) 

5U^{d) = I^Sel-aijdV (3.12) 


or, 

5U^{d) = f r^\a,,‘5e^^ + ayy5eyy + T,y5'r^ + 

JAJ—h/2 

ry,-6^y^ + T,,-S'fndV (3.13) 

Now, by using equations (2.23)-(2.27), (2.36)-(2.37),(2.57)-(2.58), and (3.4)-(3.10) in equa- 
tion (3.13), we shall get the following stiffness coefficients corresponding to the particular 
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degree of freedom as: 


= nS;) iV,(n,-) nS;) * N^{nj) 

(3.14) 

^if = ”i!!? * * ^xyirij) 

(3.15) 

Kij^ = TljJ * (^Nx{nj) * W^x "h ^xyi'^j) * + QxiP'j)) 

* "^iV d" ^xyi'IT'j) * '^,x + Qyi'^j)) 

(3.16) 

♦ Qx{nj) + * J14(%) -I- ♦ Mxy{nj) 

(3.17) 

* Qyirij) + * My{nj) + * Mxyinj) 

(3.18) 

Kif = 3.0 ♦ * Qlinj) + * M*x{nj) + * M:,{nj) 

(3.19) 

Kjp = 3.0 * nP * Qliuj) + * M*{nj) + 

(3.20) 


where from equations (2.17)-(2.22),(2.60),(3.4)-(3.10), and (3.14) -(3.20), we have 

-t-(i4i2 * WjJ + Aiz * nlj l) 

. +{Bn * + Biz * 

+ {Bi 2 * nfl + Biz * nfl ) (3-21) 

-|-(jBi4 * nf^ + Biz * nfl) 

+{Bi5*n’^^+Biz*nf}) 

-l-(Aii *w^* nfl + An * w,y * nfj) 

+{Aiz*w^* nf^ ■\-Aiz*w,y* nfl ) 

Ny{nj) = (Ai2*n$? + ^3*n2^) 

+(^22 * nfl + ^23 * nfl) 
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+(Bi 2 *ng + B23*ng) 

+(522*ng+£23*ng) (3.22) 

+(J524 * nfl + ^26 * nf^) 

+{B 25 * nf^ + B 26 ♦ n^) 

+(Ai2 *w^* nfl + A22 * w^y * nfj) 

+(A 23 * * nfl + A 23 * w,y * nf^) 

Nxy{nj) = (Ai 3 * + A 33 * 

+(A23 * nfl + A 33 ♦ nfl) 

+{Biz * + Bzz * nf^) 

+{B2z * + Biz * nf^) (3.23) 

+(Bi4 * nf^ + JB 16 * nf^) 

+(J 5 i 5 ♦ + Biz * 

+(Ai 3 * Wx * nf^ + A23 * w^y * nfj) 

+(A33 * KJ,* ♦ nfl + A 33 * w,y * nfl) 

Mx{nj) = (Bn * nf^ + Bzi * nf]^) 

+{Bi2 * Tlj J + Bzi * ^jj) 

+(An * 4 ^] + Ai 3 * rS^l) 

+( Ai 2 =*= ng + Ai 3 * 

+{Dh,*nf} + Dt,^*nf^) 

+(Ai 5 * + Aie * 

My{nj) = {Bi2*nf}+Bz2*nf]^) 

+ (^22 * A Bz 2 * ^j,*) 


Finite Element Formulation 


29 



+(f »11 • "2 + Dtj, * nj**) 
+(A,.*n£ + Dfe,.ng) 

+(A„.ng+D»„*r.S«) 

+(Cj...n£+Dfe,*n£) 

(3.25) 

^xy{^j) 

= (B,3*n<-«+B33»nS) ; 



+(^23 * + B 33 * 

+(Dts 2 * nfl + Db^ * nf^l) 

+(A34*n£+A3s*nX^) 

+ (^35 * nS + '^*>36 * ^S) 

(3.26) 





+iB 2 A*nf^+Bu*nf}) 
+iDu.*n^£ + Db,,*n^^l) 
+(A«*n£+A43*ng) 

+(I>644 * nf} + A 46 * nfl) 

+ (^45 * ^'i^v + -^646 * 

(3.27) 


= (Si5*ng + B35*n5^]) 



+(B25 * nfl + B 35 * nfl) 

+(^51 * nfl + D 653 ♦ nfl) 

+(B>652 * nfl + As3 * ^52 ) 

+ (^54 * nfl + ^656 * nf^) 

+ (^55 * 

(3.28) 
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KvK) = (Bi6*ng + 536*ng) 

■ + (B26 * nfl + B36 * ng ) 



+(B6g2 * Tlj J + Di,^ * ^j,i) 

+ (■^665 * ^^j,l 1^666 * 

(3.29) 

Qxi'^j) 

= {D,,,*nf + Ds,,*nf) 

(3.30) 


+{^sn * nf} + ^*12 * 45) 


Qy i'^j) 

1 

(3.31) 


+{Ds.^*nfl+D,,,*nf^) 


Q>i) 

= S-Ds33*nf +Z- D su*t4^) 

(3.32) 

Qyini) 

= 3 • Dsi 3 * nf -1- 3 • ♦ nf 

(3.33) 


3.1.1 Equivalent nodal load vector: 

The external work-done due to transeverse loads is expressed in terms of shape functions 

as: 

(3-34) 

i=l 

where ’e’ denotes the elemental value, and 

■jT 

= [ 0 0 5^®^ 0 0 0 0 j 


(3.35) 
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3.1.2 Non-linear equilibrium equation: 

For equilibrium of the element, we have 

= 0 (3.36) 

i.e., 

I 6 = 0 (3.37) 

Equation (3.37) gives the non-linear equilibrium equations after collecting the coeflBcients 
with displacement term and loading term separately. 


^di) = • dj - Fi (3.38) 

where, 

I Fi = Jnih^lwr} dA (3.39) 

where, I 7 is the 7x7 unit matrix. 

3.2 Solution to non-linear equilibrium equation: 

The solution algorithms for the assembled non-linear equilibrium equations (3.38) are 
based on the Newton-Raphson method, which consists of a series of linear solutions.If an 
initial estimate di, for the total displacements, gives the residual (or unbalanced) forces 
'^[di) ^ 0,then an improved value dj+i is obtained by equating to zero the linearized Taylor s 
series expansion of ^(dj+i) in the neighbourhood of dj as, 

^(dj+i) ^(dt) + Kt • ^dj = 0 (3.40) 

where Kx is known as assembled tangent stiffness matrix evaluated at d/ and is given 


by, 
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Kt — 


d<e(di) 


(3.41) 


[ m \ 

Equation (3.40)is the linear incremental equilibrium equation which gives the linearized 
approximation to the relation between the residualforces and incremental displacements, 
Adi, at a point di on the equilibrium path. The improved solution is then found as: 


dj+i = di + Adi 


(3.42) 


Equations (3. 40) and (3.42) represent the Newton-Raphson solution approach to the 
nonlinear equations (3.38). The terms '^{di) and Adi give the measures of the convergence 
of the solution. To improve the numerical stability and to give intermediate results, the load 
Fi is usually applied in steps. 


3.3 Tangent stiffness matrix: 

The tangent stiffness matrix given .in equation (3.41) may be written as. 


Kt = 




d±L 

ddi 

dd2 

•• ddm 



d^2 

ddi 

dd2 

•• ddm 



d^m 

ddi 

dd2 

•• ddm 


where m is the number of degrees of freedom per element. 
Prom equ. (3.38) and (3.41), we have 

„ dm) 


_ f.-n I I 

- gdj ° ed, 


(3.43) 


(3.44) 
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Therefore, 


KT = K,i{d) + 


f 9Kym 

\ J 


■di 


[ since |^ = 0 for conservative loads. ] 

In equation (3.45), the term 

d{Kij{d)) _d^Kij+^^Kij{d) 

m d(dj} 


(3.45) 


(3.46) 


d(Kij(d)) _ d^^Kijjd) 
d(dj) d(dj) 


(3.47) 


Therefore, from equations (3.45) and (3.47), the tangent stifiBiess matrix becomes: 


(3.48) 


We can find the non-linear part of stiffness matrix K^^{d) by synthesizing equations (3.14)- 


(3.20) and using equations (3.21)-(3.33) in these, so as: 

iVL^(i) ^ „(i)/ ^nw,Jll+Auw,^n^l 1 
1 AAizw^nfl + AiiW^nl^ J 
(lyf AziW,,nfl+Az2W,^nf^ 1 

1 +Azzw,^nf^ J 

NLrci'^) = „(2)/ \ 

1 +A23W,,nf]^ + A23W,yj^^ J 

( 2 ) f ^31 v2+^32U;,v^S 1 

1 +A33W^nfj + Azzw^^nl^ J 

{Anw^nf} + Ai3W,^n^J;l) 
+iAi2W,,n^ll + Anw,,nlJ 
+iBnW,,n^^l + Bizw^^n^^l) 


(3.49) 


(3.50) 
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+{Bi2W,^nf^ + Bi3W,,nfl) 

+{Buw,,nfl + Biew^^nf^) 

+ {Bi5W,^nf^ + 

+'^i,. • [ 

(Aiiw.ng + Ai2w,ng) 

+ Ai3W,nf^) ] 

+{A3iW,^nf} + A33W„n^^l) 

+{Bziw„nf} + B33W,^n^^l) 

+{B32W,^nf^ + B33W,^nfl) 

+{B3AW,^nf^ + B36W„ng) 

+(B35U^.,ri52 + B36W„n2) 

+^»,V • [ 

iA3iW,,nfl + A32wyf^) 

+{A33W^nf^ + A33W,nfl) ] } 

+nSJ • { (3-51) 

{A2iw„nf;l + A23W,^nf^^) 

+{A22W„nf^ + A23W,yjJ 
+{B2iW,^n^^l + B2zwy^l) 

HB22W,,nfl + B2zw,^n^^) 

+(B24W,n£+526W,n5®j) 

+(^25W„n5y^ + B2ew„nf^) 

+Wi,» ' [ 
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{A2iw,^nfl + A22W,^nfl) 
+{A2zw,^nf} + A2zw,^nf}) ] 
+{A3iw,,nfl + A33?i;„n5>^) 
+iA32W,,nfl+Az3W,^nfj 
i +iB3iW,,nfl+BzzW,,nf^) 

i 

I +{B32W, ,nf I + BzzW^^nf}) 

+{Buw,,n^^l+B3zw^,nf^) 

+{Bz5W,.n^J^ + Bz6'w,,nf}) 

+Wi,, • [ 

{Aziw,,nfl+A32w,^nfl) 

+{Az3W,^nf^ + Azzw^^nf}) ] } 

NLr.W ^ U4)| BnW,JJ!l + B2xW,^nf^ 1 

\ ABzxW^nf^ + Bzxw,^nl j 
(4>f Buw,,nfl+B2zw,yf^^ 1 

”‘”'1+B33vg+J533^„nf,, J 


NiKi^ = ^^3^.,n£+B23t/;„nS 

*’* j +B33W,.n5^^ + Bzzw,^n\^ 

(5) J Bx2W^,rifl-\-B22W,^rifl ) 

*”' \ +B32W,^nfl + B32W,^nl^ ) 


NL 1^(6) 
J^ij 


(6) f SuW^^Tlj^ + B24W,ynj^ 1 

'‘'■\+B„to,.ng + BMW,.n5„ J 

(6) I Bi«»,.n£+B26t»„"?j 1 

1 +fl36«'«nS + J 


(3.52) 


(3.53) 


(3.54) 
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NLj^m 


% 


( 7 ) 

iix 




,( 3 ) 

hx 

.(3) 


>y J>y 
3 

'i,® 


+n 


+B367^^,.n) j + ^36^, 

(7)/ ^15vS+^25Wv^i2 

+B35W„n5^^ + B35U^„nf_^ 


(3.55) 


Therefore, 


adf I “ a4^^ 

d^^Ka d^^Kii 


ad« 


84'> 


adf^ 


a4^) 


(3.56) 


and, 


dSP 


- .(1) 


= n, 




( 


+n, 


*IV 


+^33’2i!i«lJ + •^33^2 


(3.57) 


or. 




ad' 


(3) 


*>» 


+Ai3W,ng+Ai3n;,„n'_^ 


ny: + A32W 




4 - A33W 



(3.58) 


A^}_NLiAl) 

^.(3) ■ 


adv 


(3.59) 






ad)' 


{3) 


(3.60) 


a“/ff 


adf 


(3.61) 


Similarly, 
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^ , j(3) NL tAS) 

a4 - 4 

(3.62) 


^3) A - 

(3.63) 


adf ' " 

(3.64) 

and, 

dNLK§) 

a<i-p 




+’^ifv ( + ^xy(nj) • n|fj ) 

(3)/ ^ 

^'‘^*1* ’* ddf^ ’* y 

+"i.. 1, 8ij« ' S^"' - j 

(3.65) 

Therefore, 



94^ 

(3)/ N,{n,)-w,+N^y{n,)-wJ N 

1 +<’ WW ■ w,^ + N^im) ■ w, j 



ddf^ ’* Sd) ^ J 

, (3) ( dNy{ni) , , a/v„M . _ 

(3.66) 

Prom equation (3.14)-(3.20),we have 


NL tx(3) 

I\ij 

— ' '^,x + ^xyiAi) ■ '*^.v ) 



+^i?J ( • ^.v + ) 

(3.67) 


and from equation(3.21)-(3.33), we have 
dd) ' ^ 


(3.68) 
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\ ( A2inf} . nlj^^ + A22nf^ • ) 


+ « ( ^23^2 • + ^23^52 • n\^} ) 


(3.6§-) 


Thus, 


or. 


dOj 


= 5(^.nS2-nl./) + A.ng.ng) 

+ I (M2 +'433ng-n2) 




dd. 


,(3) 


1 

2 


+ x•J^, 


( 3 ) 

'i)X 


Anw,,nf} + Aiaw^n);' 
+Ai3W,,nfl + Ai3U),,n? ^ 


(3). 


w. 


+ 


+ 


1 

2 

1 

2 


+ 




nf 

*»y 


nf) 

tjj, 


+ A32W„nf,' 
+A33U^,.ng+A33W,nf_, 

,(3) i„„,„ 43) 

J>y 


( 3 ) 


ty 


>y 




ty - 


^2i^„n}j + A22t/;„nV 


+A23W,^nf^ + A-izw^^n] 

f 

+A33t^'„ni„ 


%X 

Aiiw^^nj^l + Az2W,^nf^^ 
A Azzw,^nfl 


(3.70) 


(3.71) 


ng ■ { w,. iV^(T.j) + w„ •“ Ws,("i) } ] P.72) 


where, 
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(3.73) 


(3.74) 


(3.75) 

Therefore, from equation (3.48) and (3.59), we get 

(3.76) 

Prom equations (3.48) and (3.60), we get 

(3.77) 


Prom equations (3.48) and (3.71), we get 


= K^fid) + '’‘‘Klf 

+(n£ ■ •“ «.(%) + •“ JV:^(n,)}l 

Prom equations (3.48) and (3.61), we get 

= ( 3 - 79 ) 

Prom equations (3.48) and (3.62), we get 


Prom equations (3.48) and (3.63), we get 


Kl 


= K^fid) + 


NL 1 ^( 6 ) 

Kij 


(3.81) 
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Prom equations (3.48) and (3.64), we get 

K?,] = (3.82) 


Now, by using equations (3.76)r(3.82), we can calculate the tangent stiiFness matrix 
[Kt]. 

The iterative solution procedure used in equations (3.38) and (3.40) is checked for con- 


vergence using the following criteria: 




pr.F 


1/2 


X 100 <13 


(3.83) 


where, /? is sufficiently small (~ 10 “®). 

Knowing the displacement vector {d}, we shall be able to calculate strains for the given 
laminate.Now, by using constitutive equations, we can get stresses at the certain point in 
the laminate. 


3.4 Material non-linear solution technique: 

We shall solve material non-linear problem by applying strains in increments. For the 
first value of strain^ we shall calculate the various tangent modulii, i.e.jS^ . Now we shall 
use these to update the tangent stiffness matrix [Kt] and then finding the displacement 
vector d. For these displacements, again calculating the strains through strain-displacement 
relations. Using this new set of strains we shall again calculate the new values of Et 3 ,nd 
again update the tangent stiffness matrix to get new set of displacements. This procedure 
will be repeated till the following criterion is not satisfied: 

€i+i - €i < a (3.84) 

where, € 1 * 4.1 and €{ are two consecutive values of strains during iterative procedure, and, oc is 
very small ~ Now, we shall go for the second value of the strains and shall again be 
repeating the whole procedure to find the displacement field which satisfies equation (3.84). 
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The whole procedure is repeated again for all other values of strain- vector. Thus, we 
shall be able to find the displacements upto desired accuracy. 



Chapter 4 

i 

Results and discussion 


4.1 Convergence Study 

The material properties for the different materials are shown in table 4.1.These material 
constants are taken from the [24] k [27]. For the material nonlinear problem, we have taken 
the following material properties for the material A(table 4.1) as [27]: 

=5.0 X 10®psi, Ep 22 = 2 . 0 X 10®psi-, 

A66=2.0, A44=2.0, A55=2.0, 

cTpn =0.005 X 10®psi, cTpjj =0.005 x 10®psi, 

m=2.0 


Table 4.1: Material constants of different materials. [24] k [27] 


Material 

Ell 

xlO® 

E 22 

xlO® 

Gi2 

xio® 

G 23 

xlO® 

Gi3 

xio® 

^12 

1^21 

A 

25.0 

1.0 

0.5 

0.2 

0.5 

0.25. 

0.01 

B 

30.0 

30.0 

11.4 

11.4 

11.4 

0.316 

0.316 

C 



0.3125 

0.3125 

0.3125 

0.23949 

0.24078 

D 

3.0 

1.28 

0.37 

0.37 

0.37 

0.32 

0.24 
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Table 4.2: Center Deflections of rectangular laminate using linear HSDT for different mesh 


Mesh 

Stacking 

Sequence 

Simply- 

Supported 

Clamped 

■'i- 

2x2 

o'' 

1.715807 xlO-^ 

5.83942 xlO-"* 


0^90'' 

5.644679 xlO"^ 

1.670752 xlO-^ 


0''/90''/90‘'/0° 

4.937523 xlO'^ 

1.488710 xlO-^ 


[457 - 45"L 

4.809534 xlO-i 

3.73201x10"^ 

4x4 

O'' 

1.464109 xlO-i 

5.6465 xlO"^ 


0790 " 

5.63067 xlO-^ 

1.461067 xlO"^ 

1 

079079070 " 

4.916800 xlO-^ 

1.42304 xlO-i 


[ 457 - 45 '']. 

3.292528 xlO'^ 

2.16222x10"^ 

8x8 

o'* 

1.460854 xlO-^ 

5.6167 xlO"^ 


0790 '' 

5.59695 xlO-i 

1.426537 xlO"^ 


079079070 '' 

4.889124 xlO-^ 

1.40254 xlO"^ 


[457 - 45"). 

3.265564 xlO'^ 

2.12656x10"^ 


The confidence in the finite element formulation is checked with the mesh refinements.Resiilts 
are taken for the analysis with linear HSDT for various type of laminates of material ’A’( 
Table 4.1). It is found that 4x4 mesh gives equally good results as compared to those found 
in literature. But, we have used 2x2 mesh for various calculations, as this mesh size gives 
results of admissible accuracy and less computer time, when compared with 4x4 mesh and 
8x8 mesh. Though the convergence of the values are achieved through the highly refined 
mesh, but the computational time will be too much. Therefore, we have settled with 2x2 
mesh only. The results of convergence-study for central deflection of the simply-supported 
and clamped plate of material’ A’ are shown in the table 4.2. These results show that 8x8 
mesh is quite accurate. 
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4.2 Various non-linear cases for cross ply and angle ply 
laminate 

The geometric non-linear results are taken tor the cross-ply with material (C) and for 
angle-ply laminates with material (A) and results are compared with that of [27]. The 
dimensions of the plate are a = b = 12in. and thickness h = O.OQGin. In figure 4.1, the load- 
deflection curve is drawn for 3-layer cross-ply laminate.Results are compared with those in 
[29], and, it is found that present study gives good results. The comparative study of 4-later 
cross-ply and 2-layer angle-ply laminates is shown in Figures Fig.4.2 and Fig.4.3. 

Following results can be inferred from the curves: 

i) Gross-ply laminate exhibits more non-linearity for same load and same dimensions 
when compared with angle-ply laminate. 

ii) Increasing the no. of layers reduces the effect of non-linearity. 

iii) Combined non-linear results can be best for higher loadings. 

4.3 Orthotropic plate; Load vs Deflection for a=b=12in 
and thickness h=0.138in for material(D) 

Results for the material nonlinear case are also presented for simply-supported [0"/90°/0°] 
orthotropic laminate in Fig.4.4 for a/h = 125 . Results for orthotropic plate are also taken 
for simply supported boundary conditions with various load values. Comparisons are made 
with values given in literature [24] and also with experimental values found in nonlinear 
analysis. Results are shown in Fig.4.4 and also compared with those in [24]. 

Some more results are shown in tables 4.10 and 4.11 respectively for a/h vs. central 
deflection for simply-supported, square, with material ’A’, 3-layer and 2-layer cross-ply plate. 
Results are also compared with those in [8] and [27], and further ensure the reliability of 
present formulation. 
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4.4 Cross-ply laminates: Variation of deflections with 
a/h for various HSDT models 

Results are taken for cross-ply laminates to find the variation of bending response with 
the side-to-thickness ratio, by using different combinations of non-linearities and also for 
linear HSDT analysis.lt is found that for thicker laminates any of the HSDT can be used, 

. i . 

but for thin laminates only non-linear HSDT can give admissible results.Results are shown 
for parametric study in tables 4.6, 4.7, 4.8 and 4.9 for a/h = 5, 10, 20 and 100 respectively 
and for different loadings.lt is found that by using combined nonlinear analysis the results 
are slightly more conserved than by using only material nonlinear analysis. 

Results for istropic plate are also taken for clamped boundary conditions with various 
load values. Comparisons are made with values given in literature and also with exact 
values.Results are tabulated in table 4.5. 

I 

j 

4.5 Cross-ply laminates: Deflections with polynomial 
stress-strain relations for normal stresses 

The polynomial stress-strain relations can be used in place of modified Ramberg-Osgood 
relations for the nomal stresses too. As the laminates are seldom encountered with the large 
stresses therefore we can use these polynomial equations for finding the results under the 
particular value of strains. Results are compared with those found with modified Ramberg- 
Osgood relations[27] and show equally good results. For comparison, we can see tables 4.5 


and 4.7 simultaneously. 
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Table 4.3: Material Properties of T300/5208 (pre-peg) graphite-epoxy. [36] 


Mechanical 

properties 

Values 

Strength 

properties 

Values 

Ex 

132.58GPa 

Xt 

1.515GPa 

E2 

id.SGPa 

Xc 

1.697GPa 

jE?3 

lO.SGPa 

Yt = Zt 

43.8MPa 

i Gx2 = Gxz 

5.7GPa 



' 1^12 = 1^13 

0.24 

R(shear 
strength in 
2-3 plane) 

67.6MPa 

1^23 

0.49 

S(shear 
strength in 
1-3 plane) 
=T(sh. str. 
in 1-2 plane) 

86.9MPa 


4.6 Failure analysis by using linear HSDT, geometric 
non-linear HSDT and combined non-linear HSDT. 

Results of failure analysis are found for transverse loads only.Results are taken for the 
various lamination schemes and for various combinations of non-linearities. Results are com- 
pared with those in [36]. First ply failure loads of the given laminates under the transverse 
load for the maximum strain criteria, are presented in Table(4.12), respectively for various 
theories, viz., linear, geometric nonlinear, and geometric as weU as material nonlinear com- 
bined theory. It is known that maximum strain criterion yields the lowest strength for all 
load combinations. Therefore, maximum strain criterion is used to obtain the results witm 
lower bounds. Results are compared with those found in [36]. 

In this analysis, the laminate dimensions are taken as 479.88 mm X 479.88 mm X 3.72 
mm with ply thickness as 0.2325mm([36]). The material properties axe shown in table 4.3. 
The laminate is descretized into 4X4 meshes. 
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Load p= 

^2 


Figure 4.1: Load-deflection curve for Material nonlinear case for a//i = 10. 


Table 4.4: Center Deflection for clamped, isotropic square plate for material (B),(a-b-300m, 
h=3in, Q = {qQa!^)l{Eh^) ) 


Load 

Q 

ESI 

mmm 

mm 

Present 
(linear HSDT) 

17.79 

0.237 

0.2368 

0.2417 

38.30 

0.471 

0.4699 

0.4723 

63.40 

0.695 

0.6915 

0.6961 

95.00 

0.912 

0.9029 

0.9127 

134.9 

1.121 

1.1063 

1.1217 

184.0l 

1.3231 

1.300F 

1.3243 

245.0 

1.521 

1.4928 

1.5311 

318.0 

1.714 

1.6786 

1.7160 

402.0 

1.902 

1.8555 

1.9044 





































Centre deflection, w(in.) 
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Figure 4.2: load vs.deflection for cross-ply plate. 


Table 4.5: Cross-ply: Center Deflections with polynomial stress-strain model (cubic) for nor- 



Load! 

Intensity 

20 

40 

100 

500 

. 1000 

L 

S. S.(xlO^) 

2.14996 

4.29992 

10.7498 

53.7490 

107.49920 

L 

Clamped(xlO®) 

1.26759 

2.5351 

6.337 

31.6850 

63.3700 ■ 

GNL 

S. S.(xlO®) 

2.149957 

4.29987 

10.7489 

41.2246 

69.6578 

GNL 

Clamped(xlO^) 

1.267578 

2.53504 

6.3356 

26.4623 

42.6645 

MNL 

s. s.(xi(y') 

2.14802 

4.298014 

10.744 

44.5612 

75.8745 

MNL 

Clamped(xlO®) 

1.26712 

2.53483 

6.33568 

28.5585 

44.6792 

GMNL 

S. S.(xlO^) 1 

2.147864 

4.297966 

10.7^^ 

43.4532 

72.6547 

GMNL 

Clamped(xlO^) i 

1.26687 ! 

2.52475 

6.3153 " 

27.9735 

i • ^ 1 

43.2375 
^ 1- 


GMNL=geometric and material, combined nonlinear.) 









































Nondimensionalized center def. , w — w / h 


j^.6 Failure analysis by using linear HSDT, geometric non-linear HSDT and combined non-linear 
HSDT. 49 


Clamped Angle Ply (45 / -45 ) 



Load parameter, P 


Figure 4.3: load vs.defiection for angle ply laminate. 



Load 

Intensity 





1000 

L 

S. S.(xlO^) 

3.034289 

6.06857 

15.17144 

75.85720 

151.71450 

L 

ClamDedfxloT 

2.3310 

4.6621 

11.1655 

55.8275 

111.65520 

GNL 


3.0342886 

6.068573 

15.17136 

56.4355 

92.8843 

GNL 


2.33107^ 

4.66213^ 

11.165501 

47.9756 

78.5465 

MNL 


3.034288 

6.068575 

15.17139 

59.7865 

95.2374 

MNL 


2.331073 

4.662121 

11.16548 

49.5575 

81.6534 

GMNT. 


3.034288 

6.068574 

15.171375 

57.6785 

94.2134 


Clamped(xlO^ 

2.331076 

4.662142 

11.61519, 

48.5764 

79.8745 


/ — iliiuax, 

GMNL=geometric and material, combined nonlinear.) 
















































Centre deflection, w (in.) 
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O 


Figure 4.4: load vs.deflection for orthotropic plate. 


Table 4.7: For afh = 10 



Load 

Intensity 

20 

40 

100 

500 

1000 

L 

S. S.(xlO^) 

2.14996 

4.29992 

10.7498 

53.7490 

107.49820 

L 

Claniped(xlO^) 

1.26759 

2.5351 

6.337 

31.6855 

63.3720 

GNL 

S. S.(xlO^) 

2.149957 

4.29987 

10.7489 

44.8372 

66.4735 

GNL 

Clamped(xlO^) 

1.267578 

2.53504 

6.3356 

23.1563 

39.8745 

MNL 

S. S.(xlO^) 

2.14992 

4.299612 

10.744 

47.5637 

68.76540 

MNL 

Clamped(xlO^) 

1.26758 

2.53506 

6.33598 

26.7635 

43.5675 

GMNL 

S. S.(xlO^) 

2.149972 

4.299998 

10.75093 

45.7865 

67.9725 

GMNL 

Clamped(xlO^) 

1.26757 

2.53503 

6.3355 

24.8747 
• *• 

41.9736 


(where, L=linear, GNL=geometric non-linear, MNL=niaterial non-linear. 


GMNL=geometric and material,combined nonlinear.) 




4-6 Failure analysis by using linear HSDT, geometric non-linear HSDT and combined non-linear 
HSDT. e 



Figure 4.5: Load-deflection curve for Material nonlinear case for a/h = 5. 


Table 4.8: For a/h = 20 



Load 

Intensity 

20 

40 

100 

500 

1000 

L 

S. S. (xlO^) 

2.4316 

4.8632 

12.1580 

60.7902 

121.58020 

L 

Clamped (xlO^) 

8.2784 

16.5569 

41.3924 

206.9620 

413.9240 

GNL 

S. S. (xlO^) 

2.43147 

4.86209 

10.1403 

45.3213 

70.8735 

GNL 

Clamped (xlO^) 

8.2771 

16.5464 

38.2288 

168.9427 

288.3465 

MNL 

S. S. (xlO^) 

2.431469 

4.86205 

10.1399 

49.8745 

73.5467 

MNL 

Clamped (xlO^) 

8.27746 

16.5487 

38.2656 

173.9473 

294.3475 

GMNL 

S. S. (xlO^) 

2.431473 

4.862093 

10.14032 

47.1392 

71.9685 

GMNL 

Clamped (xlO^) 

8.27703 

16.5452 

38.2113 

171.7438 

291.7385 


(where, L=linear, GNL=geometric non-linear, MNL=material non-linear, 


GMNL=geometric and material,combmed nonlinear.) 




LiRRARt 
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Table 4.9: For o/A = 100 



Load 

Intensity 

20 

40 

100 

500 

. 1000 

L 

■ S. S. 

13.41730 

26.83461 

67.086353 

335.431765 

670.86353 

L 

Clamped 

3.20981 

6.41962 

16.04905 

80.24525 

160.49050 

GNL 

S. S. 

4.706729 

6.317180 ■ 

8.994607 

8.999920 

8.999945 

GNL 

Clamped 

1.23547056 

1.65246093 

2.35682114 

2.368799 

2.37 23522 

MNL 

S. S. 

4.78049241 

6.45674941 

9.2899640 

9.355345 

9.3702347 

MNL 

Clamped 

1.273988 

1.745240 

2.6299976 

2.6498564 

2.661042 

GMNL 

S. S. 

4.67617653 

6.26272171 

8.8931193 

9.116382 

9.286639 

GMNL 

Clamped 

1.21966965 

1.618488 ^ 

2.276312 

2.5693754 

2.628495 


(where, L=linear, GNL=geometric non-linear, MNL=material non-linear, 
GMNL=geometric and material,conibined nonlinear.) 


Table 4.10: Nondimensionalised center deflections for (0°/90®/0‘’) cross-ply S.S. square plate 
under UDL {qq = 1.0;w = '“^fJ^- ]Material{A)) 


afh 

w ([8]) 

w(present) 
(Geometric- 
nonlinear HSDT) 

2 

7.7661 

7.7660 

4 

2.9091 

2.90907 

10 

1.0900 1 

1.0900 

20 

0.7760 

0.77594 

50 

0.6838 

0.68376 

100 

0.6705 

0.67046 
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Center deflection for (0°/90°) cross-ply S.S. square plate under UDL (go = 1.0;ty = 
Material = (^)) 

Table 4.11: Deflection for simply-supp. boundary conditions of square plate material (A). 


i 


afh 

w([271) 

w(present) 
(linear HSDT) 

5 

2.5813 

2.5802 

EM 

1.9192 

1.9184 

Jim 

■fTOIr/il 

1.7398 

100 

WTiW^l 

1.6577 


Table 4.12: First ply failure results for S.S. case.[36] 


Stacking sequence 

failure 
load for 
linear 

HSDT 

failure 
load for 
geometric 
nonlinear 
HSDT 

failure 
load for 
combined 
nonlinear 
HSDT 

Max. 
deflection 
for combined 
nonlinear 
HSDT 

Kfgn 

Failed 
guass pt. 
no. 

(+45/ _ 45/0/90)2. 

1610 

1569 

1547 

60 

16 

144 

(■+•45/ - 45 / 0 / 0 ) 2 . 

1530 

1478 

1455 

3.72 

16 

144 

(+45/ - 45 ) 4 . 

1470 

1440 

1425 

4.18 

16 

144 

( 0 / 90 ) 4 . 

2066 

2040 

2010 

4.32 

16 

106 










































Chapter 5 

i 

i 

Summary and Conclusions 


5.1 Summary 

The aim of present investigation has been to study the nonlinear response and the first 
ply failure of laminated plates under transverse loading. The overall investigation is based on 
the finite element formulation using the higher order shear deformation theory with a nine 
noded Lagrangian element having seven degrees of freedom per node. Green’s strain vector 
imposed by Von Karman’s assumptions is used to incorporate geometric nonlinearity.The 
governing finite element equations have been derived by using the total potential energy 
approach and then assembled to obtain a set of nonlinear algebraic equations.The equations 
have been solved by using Newton-Raphson method using incremental cum iterative pro- 
cedure.The material non-linearity is also incorporated by using modified Ramberg-Osgood 
relations for normal stresses and for shear strains a cubical formula is used.The tangent 
modulii for various directions axe obtained for each strain level and then stiffness matrix is 
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updated for each strain level. Picard’s direct iterative procedure is used to solve the material 
non-linear case. For geometric and material nonlinear combined case, both these iterative 
techniques are used.Different types of laminates have been analysed for the simply-supported 
boundary conditions for maximum strain criterion only.The first failed layer number, the first 
failed guass point no. and the associated transverse deflections have also been found. 

5 . 2 Conclusions 

On the basis of the present study, the following important conclusions are made: 

• The cross-ply laminates show more non-linearity than angle-ply laminates. 

• By using combined non-linearity, cross-ply laminates behave to be stiffer than angle-ply 
laminates but also with increased non-linearity. 

• The polynomial stress-strain relations can be used for the normal stress-strain case also 
instead of using the modified Ramberg-Osgood equations under the limits of strains 
not more than a particular value.Thus a good quantity of computational time can be 
saved. 

• The cross-ply laminates has the largest strength among the various laminates consid- 
ered. 

• In general,the failure occurs near the corner of the laminate. 

• The maximum transverse deflection under pure transverse load is greater for the linear 
HSDT than the geometric nonlinear HSDT. 

• The maximum transverse deflection under pure transverse load is lesser for the com- 
bined non-linear case. 
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5.3 Scope for future work 

Some problems of interest are: 

• Buckling analysis with combined non-linearity. 

• Dynamic analysis of laminates with combined loading. 

• Dynamic stability analysis for combined non-linearity. 

• Various analyses for different types of loadings. 

• Various analyses under thermal loads with combined non-linearity. 

• Progressive failure analysis of laminates under combined loads for combined non- 
linearity. 

• Effect of presence of cut-outs on the strength and failure of laminates with combined 
non-linearity. 
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